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ABSTRACT 



1. INTRODUCTION 



The nature of most of the ~ 300 high-energy 7-ray 
sources discovered by the EGRET instrument aboard 
the Gamma-ray Observatory (GRO) between 1991 and 
1999 is one of the greatest enigmas in high-energy as- 
trophysics. While about half of the extragalactic sources 
have been optically identified with Active Galactic Nuclei 
(AGN), only a meagre 10% of the galactic sources have 
a reliable identification. Interestingly, a few of them (7 
in total) have been identified with Isolated Neutron Stars 
(INSs) thanks to the coincidence with known pulsars and 
to the discovery of 7-ray pulsations at the expected pe- 
riod. The low success rate in the identification of 7-ray 
sources in the crowded regions of the galactic plane has 
mainly to be ascribed to the local crowding of potential 
optical counterparts and to the large 7-ray error boxes 
(of the order of one degree in radius) which prevented 
a straightforward optical identification. Indeed, a multi- 
wavelength identification strategy, based on a systematic 
coverage of the 7-ray error boxes, has been the only do- 
able approach. 

The situation is now greatly improving thanks to the ob- 
servations performed by the Fermi Gamma-ray Space 
Telescope (launched in June 2008) which, thanks to the 
LAT instrument, provides a factor of 50 improvement in 
sensitivity and a factor of 10 improvement in positional 
accuracy. However, while the smaller error boxes will 
make the multi-wavelength follow-ups easier, the larger 
sensitivity will enormously increase the number of de- 
tected 7-ray sources, requiring an even larger effort in 
the multi-wavelength follow-ups. This effort can not be 
obviously sustained by targeted observations only and it 
would greatly benefit from multi-wavelength data and ad- 
vanced data products available world wide through the 
science data centres and interfaced by the Virtual Ob- 
servatory (VO) tools. In this contribution, I outline the 
science case, the multi-wavelength observation synergies, 
and the requirements for both the the science data centres 
and the VO. 
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Launched in 1991, the Gamma-ray Observatory (GRO) 
was the second of the NASA's Great Observatories af- 
ter the Hubble Space Telescope (HST). In less than eight 
years of operational lifetime, the GRO performed the 
deepest mapping ever of the 7-ray sky at energies above 
50 MeV with the EGRET instrument. In particular, it dis- 
covered almost 300 new 7-ray sources (Hartman et al. 
1999), more or less equally distributed between extra- 
galactic > 10°) and galactic (|6| < 10°), correspond- 
ing to almost a factor of 100 improvement with respect to 
the results of the previous high-energy 7-ray observatory, 
the ESA's COS-B satellite. GRO observations clarified 
the nature of the sources detected at high-latitude, half of 
which have been identified with different classes of Ac- 
tive Galactic Nuclei (AGN), thus joining the quasar 3C 
273, the first extragalactic 7-ray source discovered and 
identified by COS-B. However, only 10% of the galactic 
7-ray sources were identified during and after the GRO 
mission (see Caraveo 2008, for a recent review). This 
meagre success rate is mainly to be ascribed to the large 
error box of the 7-ray sources (of the order of one de- 
gree in radius) which prevented the straightforward op- 
tical identification in the crowded regions of the galactic 
plane. 

The nature of the unidentified galactic 7-ray sources thus 
remained a mystery after the end of the GRO mission 
and it has been subject to many speculations. Interest- 
ingly, most of the few identified galactic sources turned 
out to be Isolated Neutron Stars (INSs), unambiguously 
identified as such by the discovery of 7-ray pulsations 
following the positional coincidence with a radio pul- 
sar. It is thus natural to assume that many other uniden- 
tified galactic 7-ray sources might be INSs. However, 
it most cases the poor 7-ray statistics and the lack of a 
positional coincidence with a known radio pulsar ham- 
pered the search for a periodicity which would automati- 
cally certify these sources as INSs. Indeed, some INSs (if 
not the majority) may actually turn out to be radio-silent, 
like Geminga, the first INS of this class to be discovered 
(Bignami & Caraveo 1996 and references therein), which 
would obviously depreave the astronomer from a radio 
reference period for the timing analysis. Since the large 
7-ray positional uncertain heavily affects the reconstruc- 
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tion of a coherent pulsed signal, this is almost certainly 
unrecognized in a "blind" periodicity search, which ham- 
pered many identification efforts in the past years. Of 
course, all known classes of galactic X-ray emitters, e.g. 
X-ray binaries, microquasars, magnetars, supernova rem- 
nants (SNRs), are also viable counterparts to some 7-ray 
sources. In addition, it can not be excluded that some 7- 
ray sources might be indeed associated with entirely new 
classes of astrophysical objects, whose existence passed 
unnoticed so far. Apart from INSs which can be identified 
via the detection of 7-ray pulsations, by itself not always 
feasible, no other class of galactic high-energy sources 
features an evident identification signature. Thus, in most 
cases 7-ray source identification studies can be pursued 
only through a systematic multi-wavelength approach. 



2. THE MULTI-WAVELENGTH IDENTIFICA- 
TION APPROACH 



The multi-wavelength identification approach, success- 
fully tested in the case of Geminga (see Bignami & Car- 
aveo 1996 for a summary), passes through a sequence of 
"top-down" steps, i.e. consisting of observations at in- 
creasing wavelengths and angular resolution. The first 
step is the systematic mapping of the 7-ray error boxes 
through X-ray observations to pinpoint possible X-ray 
counterparts to the 7-ray source, under the assumption 
that a 7-ray source should also be detectable at its "near- 
est neighbour" wavelengths. The second step is the map- 
ping of the detected X-ray source positions with opti- 
cal/IR observations to search for potential counterparts. 
This yields to the X-ray source identification and classi- 
fication and thus to the selection of the best X-ray source 
candidate counterpart to the 7-ray source, which ulti- 
mately yields to its identification. 

The X-ray source identification process is complex and 
requires different types of information. From the X-ray 
side, classification evidence comes from the X-ray spec- 
trum, from the derived hydrogen column density Nh 
which, together with the X-ray coordinates, gives an in- 
dication on the source distance and thus on whether it 
is galactic or extragalactic, from the X-ray flux hardness 
ratio (HR), which is characteristic of the X-ray source 
class, and, of course, from long/short-term X-ray vari- 
ability and periodicity. Another parameter typical of the 
X-ray source class is the ratio between the X-ray and op- 
tical flux F x /F opt (see, e.g. Voges et al. 1999). On 
the optical side, classification evidence comes from the 
spectral energy distribution (SED) of the optical counter- 
parts, which requires either long slit spectroscopy (for the 
brightest objects) or a photometry coverage in at least 5 
passbands for a reliable classification, from optical vari- 
ability, from an optical proper motion, which obviously 
suggests a galactic source, and from the object morphol- 
ogy (stellar or extended). 

Of course, barring the large amount of observing time 
involved to map a full 7-ray error box in both X-rays 
(~ 150 ks) and in the optical (~ 80 hours for five pass- 
bands), the identification of a hundred (or more) X-rays 
sources by assembling the collected X-ray/optical infor- 



mation through a human supervised decision tree is ob- 
viously a time consuming process. Because of that, the 
multi-wavelength identification approach has been ap- 
plied only for a few, well-selected, cases like, e.g., in La 
Palombara et al. (2006) where coordinated X-ray/optical 
observations with XMM-Newton and the ESO/MPG Wide 
Field Imager (WFI) allowed the identification of likely 
INS counterparts for the 7-ray sources 3EG J0616— 3310 
and J1249— 8330, capitalizing on their better positional 
accuracy (« 0.5 degrees) and on the 33' x 33' WFI field 
of view which well matches that of the EPIC instrument 
aboard XMM-Newton. Similar coordinated X-ray /optical 
observations yielded to the identification of a handful 
more 7-ray sources (see La Palombara et al., 2006, and 
references therein). As it can be expected, a systematic 
multi-wavelength follow-up for, e.g. all the ~ 80 uniden- 
tified galactic GRO 7-ray sources with error box smaller 
than one degree has never been tried since it would have 
required at least 10 Ms of XMM-Newton time to obtain 
a reasonably deep X-ray coverage of all error boxes, and 
the equivalent of one year of observations with a 2.5m- 
class survey telescope. Such an amount of observing 
time, pipeline data reduction load, and data interpretation 
overheads is only affordable to large, world-wide collab- 
orations. 



3. THE AGILE AND FERMI GAMMA-RAY MIS- 
SIONS 



Perspectives for 7-ray astronomy became brighter with 
the launch of the AGILE and GLAST, now Fermi Gamma- 
ray Space Telescope, satellites. AGILE, launched on 
April 2007, launched by the Italian Space Agency (ASI) 
mounts a main instruments for high-energy 7-rays: the 
Gamma-ray Imaging Detector (GRID) which covers the 
30 MeV-50 GeV spectral range with a sensitivity equal 
(or slightly better) than that of the GRO/EGRET. The 
GRID instrument has a field of view of 3 steradiants 
(about 1/4 of the sky) and a positional accuracy better 
than 0.3 degrees, far better than that of EGRET. Further- 
more, AGILE is equipped with the Super-AGILE detec- 
tor which works in the hard X-ray range (10-40 keV), 
with a smaller field of view (0.8 steradiants) but a posi- 
tional accuracy better than 3', which makes it extremely 
useful for prompt or parallel coverage of the 7-ray error 
boxes. For a detailed description of the AGILE mission 
see Tavani et al. (2008). FermQ, launched by the NASA 
on June 2008, mounts the Large Area Telescope (LAT) 
which has a factor of 50 higher sensitivity with respect 
to the GRO/EGRET. Compared to the AGILE/GRID, the 
LAT instrument has a larger spectral coverage (20 MeV- 
300 GeV), a comparable field of view (2.5 steradiants), 
but a positional accuracy better than 0.1 degrees. In the 
first year of operation Fermi will perform a deep all-sky 
scan (3 hours per single scan) which should lead to the de- 
tection of about 10000 7-ray sources with a significance 
better than 5a. 

Among the early results, AGILE observations yielded to 

1 http://fermi.gsfc.nasa.gov/ 
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the identification of the GRO source 3 EG J2021+3716 
with the radio pulsar PSR J2021+3651 (Halpern et al. 
2008). More spectacularly, Fermi has identified a new 7- 
ray pulsar (Abdo et al. 2008) with the unidentified source 
3 EG J00 10+7309 in the CTA 1 SNR. This was possible 
only thanks to the unprecedented positional accuracy of 
the LAT instrument, crucial to reduce pulsar timing prob- 
lems, which has allowed to detect pulsations from INSs 
even without the aid of a reference radio period. At the 
time of writing this contribution, Fermi has already dis- 
covered that 1 1 more, so far unidentified, 7-ray sources 
are indeed 7-ray pulsars which, like that in the CTA 1 
SNR, are radio-silent. This is a further evidence that 
many unidentified galactic 7-ray sources are likely radio- 
silent INSs, as put forward a long ago since the identifica- 
tion of Geminga in the early 1990s. However, it has to be 
reminded that while many unidentified 7-ray sources can 
be straightly identified as INSs by Fermi through the de- 
tection of 7-ray pulsations, some of them might be either 
too faint for a periodicity search or might not be INSs but 
other classes of astronomical objects for which no timing 
signature is expected. Thus, multi-wavelength observa- 
tions will still play a crucial role in the identification of 
the newly detected 7-ray sources. 



4. EXPLOITING THE MULTI-WAVELENGTH 
ARCHIVES AND THE VO 



The multi-wavelength identification of serendipitous 7- 
ray sources detected by Fermi in the first year all-sky scan 
is a major project of the LAT collaboration which is co- 
ordinating dedicated follow-ups both in the X-rays, e.g. 
with XMM and Swift, and in the optical, e.g. with the 
Gemini and the VLT. While the improved positional ac- 
curacy of the LAT instrument (<0.1 degrees) results in 
much smaller fields to be covered by multi-wavelength 
mapping, thus reducing the follow-up observation load 
for a given source, its unprecedented sensitivity results in 
many more detected sources, thus increasing the number 
of fields to be mapped. This means that any large scale 
7-ray source identification campaign will require an even 
more huge multi-wavelength follow-up work which can 
not be accomplished only through a programme of tar- 
geted observations. 

In this respect, an important resource is represented by 
the huge amount of data available in the multi-wavelength 
archives maintained by world-wide science data centres, 
which help to reduce the observation loads. However, the 
level of the data available in different archives is not even. 
In the X-rays, post operational archives (POAs) are avail- 
able for the ROSAT and ASCA missions and incremental 
archives are available for Chandra, XMM-Newton, and 
Swift, in all cases providing the user with fully processed 
and calibrated data sets. In addition, in several cases 
archival X-ray data sets come with associated object cat- 
alogues with basic information on the X-ray source prop- 
erties (position, flux, spectrum, extension, etc.). This is 
the case of the ROSAT catalogue^] of the all-sky survey 



and of the pointed observations. Recently, X-ray cata- 
logues have been produced both for Chandra, an early 
release based on the first three years of operations (Ro- 
mano et al. 2008; Romano et al., these proceedings), and 
for XMM-Newton (Watson et al. 2009). Indeed, the lat- 
ter (a.k.a. 2 XMM), is the largest X-ray catalogue ever 
produced, with a total of k, 200000 unique sources (i.e., 
accounting for multiple detections) over more than 500 
square degrees. At the same time, the compilation of the 
first Swift X-ray source catalogue is under way. In the fu- 
ture, this already huge database will be increased by ob- 
servations performed by the eROSLTA X-ray telescopeQ, 
to fly in 201 1 as part of the payload of the french/russian 
Spectrum-X-Gamma satellite. With a field of view of 
1° x 1° and an improved sensitivity to higher energies 
eROSLTA will perform the first X-ray all sky survey in 
the 0.1-12 keV band, yielding at least a factor of 10 in- 
crease in limiting flux with respect to the ROSAT all sky 
survey. While the wealth of databases of X-ray sources is 
certainly a major advantage, the vast majority of them is 
still unclassified, which represents the major show stop- 
per in the multi-wavelength 7-ray source identification 
flow. There are several reasons for that. 
The first reason is that accurate X-ray source identifica- 
tion surveys, based on optical spectroscopy follow-ups, 
are affordable for a few selected fields only and are lim- 
ited to the brightest field objects. On the other hand, iden- 
tification surveys on larger scales, and down to fainter 
flux limits, can be performed based on multi-band imag- 
ing photometry and the associated object catalogues. In 
this case, VO tools would be important to perform more 
customised object matching between X-ray and optical 
catalogues using, e.g. different matching radii depend- 
ing on the quoted positional accuracy of the X-ray and 
optical sources. This would decrease the number of mis- 
matches and it would easy the identification process. It is 
clear that an identification approach based on catalogue 
matching can be handled only if automatic X-ray source 
classification tools are available. Assembling such tools, 
based on a human unsupervised decision tree algorithm 
is obviously not easy. A very interesting prototype tool 
for X-ray source classification was ClassX (McGlynn et 
al. 2002) which, however, was originally fine-tuned for 
ROSAT sources only. A more advanced tool has clearly 
to be versatile enough to adapt to X-ray sources detected 
by different satellites and instruments since some clas- 
sification parameters, like the computed X-ray flux and 
the hardness ratio, critically depend on the assumed en- 
ergy band. Future X-ray source classification tools could 
take advantage of already available optical classification 
tools developed by the VO scientists, like e.g. VOSpec 
(Baines et al. these proceedings) which classifies an ob- 
ject based on its observed SED and the comparison with 
model spectral libraries. Further advancements for such 
tools would be to allow SED fits to account for variable 
interstellar extinction and to account for different colour 
transformations between the photometry of different ob- 
ject catalogues, derived from observations taken with dif- 
ferent instrument and filters, calibrated with respect to 
different photometric systems, and sometimes expressed 



2 ttp://www.mpe. mpg.de/xray/wave/rosat/catalogue/index.php 



3 http://www.mpe.mpg.de/projects.html#erosita 
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in different units (e.g. Vega magnitudes or AB magni- 
tudes). This is crucial if, for instance, one would like 
to derive a consistent SED from the photometry of the 
SDSS, the GSC-2, and 2MASS catalogues. 
The second reason is that most of the currently available 
public optical databases do not provide an adequate sup- 
port for a systematic X-ray source identification work. 
Indeed, so far most of these works have been carried 
out using object lists matched from the GSC-2, USNO- 
B1.0, and the SDSS in the optical, which have a lim- 
iting magnitude of B ~ 22, and from 2MASS in the 
near-infrared (NIR), which has a limiting magnitude of 
Ks ~ 15. Of course, this represents a severe limitation 
since the deeper flux limits now reached by the available 
X-ray catalogues requires similarly deeper optical cata- 
logues to homogeneously sample the Fx / F opt parame- 
ter space. Thus, much deeper optical/NIR object cata- 
logues are needed. An important step forward will be 
made thanks to the new optical/NIR sky survey which 
will be performed with dedicated 4m-class survey tele- 
scopes like, in the southern hemisphere, the ESO VST 
and the ESO/UK VISTA. Survey data will be fully re- 
duced/calibrated by data reduction pipelines and multi- 
band object catalogues will be produced, ready to be 
used. For the XMM-Newton sources, a valid support will 
come for the recently released serendipitous catalogue 
of objects detected by its Optical Monitor (OM) tele- 
scope, which provide multi-band optical and near ultra- 
violet (NUV) photometry down to B ~ 24.5 for the 
longest pointings (see Still et al. these proceedings). Of 
course, much deeper serendipitous surveys can be built 
from observations performed with 8m-class telescopes, 
like the VLT, the Gemini, and with the HST. Although 
such surveys would be mostly incomplete in terms of 
colour coverage, they would certainly represent a valu- 
able support for at least some X-ray identification works. 
Unfortunately, only rarely the science data centres pro- 
vide advanced data products, e.g. fully reduced and cali- 
brated frames, image stacks, and object catalogues which 
would be needed for this project. This is the case, for 
instance, for the HST for which science-ready, on-the-fly 
re-calibrated, imaging data and image stacks are available 
through the STScI, ST-ECF, and CADC web sites. Next 
to come it is probably a multi-instrument HST object cat- 
alogue. The availability of advanced data products would 
certainly be an important commitment for other science 
data centres. For, e.g. the VLT the possibilities of ac- 
cessing advanced data products through the ESO science 
archive are in perspective very good since each imag- 
ing instrument is supported by dedicated data reduction 
pipelines and, in most cases, by object detection tools. 



5. SUMMARY 



So far, only about half of the ~ 300 7-ray sources dis- 
covered by the NASA GRO satellite (1991-1999) have 
been identified, and the rate reduced to less than 10% for 
the galactic ones. Thus, the nature of the galactic 7-ray 
sources has been so far one of the greatest enigmas in 
high-energy astrophysics. Solving this enigma is one of 



the major science goals in the next decade. This goal will 
be pursued starting from observations performed by the 
new generation of 7-ray observatories: AGILE (launched 
in 2007) and especially Fermi (launched in 2008). Apart 
from a better characterisation of sources already detected, 
many of which will be straightly identified thanks to the 
improved statistics, Fermi observations will yield to the 
detection of thousands more 7-ray sources which will 
await for identification. In many cases, this will be 
feasible only through a coordinated, systematic multi- 
wavelength approach. On a large scale, this is a huge task 
which is unsustainable for targeted observations. How- 
ever, it will certainly benefit of the enormous amount of 
data available in world wide archives. Of course, a bet- 
ter scientific exploitation of this data requires the help of 
both the various science data centres to provide the much 
needed advanced data products, like reduced/calibrated 
data and associated object catalogues, and of the VO to 
provide data analysis tools to make the matching of multi- 
wavelength catalogues easier and to provide automated 
object classification tools. 
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